SUB-FINSLER GEOMETRY IN DIMENSION THREE 
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Abstract. We define the notion of sub-Finsler geometry as a natural generalization of sub- 
'^i.' ' Riemannian geometry with applications to optimal control theory. We compute a complete 

f^ , set of local invariants, geodesic equations, and the Jacobi operator for the three-dimensional 

^^ ■ case and investigate homogeneous examples. 

1. Introduction 

CN , Much attention has been given in recent years to sub- Riemannian geometry; it is a rich 

subject with many apphcations. In this paper we introduce the notion of sub-Finsler geom- 

^ ■ etry, a natural generahzation of sub-Riemannian geometry. 

Q . The motivation for this generahzation comes from optimal control theory. A control 

(-H ! system is usually presented in local coordinates as an underdetermined system of ordinary 

differential equations 

(1.1) x = f{x,u), 

where x G M" represents the state of the system and m G M'' represents the controls, i.e., 
►^ . variables which may be specified freely in order to "steer" the system in a desired direction. 

Q\ [ More generally, x and u may take values in an n-dimensional manifold X and an s-dimensional 

£^ ' manifold U, respectively. Typically there are constraints on how the system may be "steered" 

\^ ■ from one state to another, so that s < n. The systems of greatest interest are controllable, 

O . i.e., given any two states Xi,X2, there exists a solution curve of (|1.1|) connecting xi to X2. 
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Consider the large class of systems which are linear (but not afiine linear) in the control 
variables u and depend smoothly on the state variables x, i.e., systems of the form 

(1.2) X = f{x)u, 

B : where fix) is a matrix whose entries are arbitrary smooth functions of x. This class is by no 

^ ! means all-inclusive, but it does contain many systems of interest; an example is given below. 

'k>( I For such a system, admissible paths in the state space are those for which the tangent vector 

^ • to the path at each point x G X is contained in the subspace D^ C T^X determined by the 

image of the n x s matrix f{x). Often this matrix is smooth and has constant rank s, in 
which case D is a rank s distribution on X. (In this case the variables (x, u) may be regarded 
as local coordinates on the distribution (X, D).) Thus the admissible paths in the state space 
are precisely the horizontal curves of the distribution D, i.e., curves whose tangent vectors 
at each point are contained in D. By a theorem of Chow [7j, the system (jl.2p is controllable 
if and only if the distribution D on X is bracket-generating, i.e., if the iterated brackets of 
vector fields contained in D span the entire tangent space at each point x G X. 
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Given a distribution (X, D) representing a system of the form ()1.2j) . we next consider the 
problem of optimal control: what is the most efficient path between two given points in X? 
In order to answer this question, we must have some measure of the cost required to move 
in the state space. This measure is typically specified by a first-order Lagrangian functional 
L defined on the horizontal curves of D: given a horizontal curve 7 : [a, b] -^ X, the action 
£(7) is defined to be 

-^(7)= / L{x,x) dx = / L{x,u)dx, 

J ^ J 'y 

where, since 7 is a solution curve of ()1.2|) . we define L{x,u) = L{x, f{x)u). Often the 
Lagrangian has the form 

L{x,u) = J gij{x)u'^u^ 

(summation on repeated indices being understood), and in this case it defines a sub-Rieman- 
nian metric {,) on D (i.e., a Riemannian metric on each subspace D^ C T^X) in the obvious 
way. Horizontal paths which minimize the action functional are precisely the geodesies of 
the sub- Riemannian metric. 

Example 1.1. Consider a wheel rolling without slipping on the Euclidean plane E^. The 
wheel's configuration can be represented by the vector *(x, y, (p, ip), where (x, y) is the wheel's 
point of contact with the plane, is the angle of rotation of a marked point on the wheel 
from the vertical, and ijj is the wheel's heading angle, i.e., the angle made by the tangent 
line to the curve traced by the wheel on the plane with the x-axis. Thus the state space has 
dimension four and is naturally isomorphic to M."^ x S^ x S^. 

The condition that the wheel rolls without slipping is equivalent to the statement that its 
path ^{x{t),y(t),(p(t),ip{t)) in the state space satisfies the differential equation 
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for some control functions Ui{t), U2{t). Thus the velocity vector *(i;, y, 0, ip) of any solution 
curve must lie in the distribution D spanned by the vector fields 

A natural sub- Riemannian metric on D is obtained by declaring the vector fields Vi, V2 to 
be orthonormal, i.e., by setting 

{uiVi + U2V2, uiVi + U2V2) = ul + u\. 

The integral of this quadratic form measures the work done in rotating the heading angle if) 
at the rate ip and propelling the wheel forward at the rate ip. 

But what if the natural measure on horizontal curves is not the square root of a quadratic 
form? For instance, suppose we modified Example 11.11 by rolling the wheel on an inclined 
plane? (Assume that the wheel has sufficient friction to remain motionless if no energy 
is put into the system.) We would expect more energy to be required to move the wheel 
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uphill than downhill, so the natural Lagrangian would not even be symmetric in u (i.e., 
it would not satisfy the condition L{x,—u) = L{x,u)), let alone be the square root of a 
quadratic form in u. It is not difficult to imagine examples where the dependence of L on m 
becomes quite complicated as u changes direction. This leads us to generalize the notion of 
a sub-Riemannian metric on (X, D) by replacing the Riemannian metric on each subspace 
Dj; C TxX with a Finsler metric. 

Recall that a Finsler metric on a manifold M is a function 

F : TM ^ [0, oo) 

with the following properties: 

(1) Regularity: F is C°° on the slit tangent bundle TM \ 0. 

(2) Positive homogeneity: F{x,Xy) = XF{x,y) for all A > 0. (Here x is any system of 
local coordinates on M and {x, y) is the corresponding canonical coordinate system 
on TM.) 

(3) Strong convexity: The n x n Hessian matrix 



dy'' dy^ 

is positive definite at every point of TM \ 0. 
(For details, see jl].) In other words, a Finsler metric on a manifold M is a smoothly varying 
Minkowski norm on each tangent space Tj-M. 

Condition El implies that the "unit sphere" in each tangent space T^M (also known as 
the indicatrix for the Finsler metric on T^M.) is a smooth, strictly convex hypersurface 
enclosing the origin 0^. G T^^M. The converse is almost - but not quite - true: there exist 
strictly convex hypersurfaces for which the corresponding Hessian matrix is only positive 
semi-definite along a closed subset; see P for examples. We will say that a hypersurface 
Sj. C Tj-M which encloses the origin is strongly convex if it is the indicatrix for a Minkowski 
norm on T^M; thus strong convexity implies strict convexity, but not vice-versa. 

In the Riemannian case, the indicatrix must be an ellipsoid centered at 0^:, but in the 
Finsler case it may be much more general. In particular, it need not be symmetric about 
the origin. 

We are now ready to define our primary object of study. 

Definition 1.2. A sub-Finsler metric on a smooth distribution D of rank s on an n- 
dimensional manifold X is a smoothly varying Finsler metric on each subspace D^ C T^X. A 
sub-Finsler manifold, denoted by the triple (X, -D, F), is a smooth n-dimensional manifold X 
equipped with a sub-Finsler metric F on a bracket- generating distribution D of rank s > 0. 
The length of a horizontal curve 7 : [a, 6] — *> X is 

£(7)= / F(7(t))rft. 



Replacing the Riemannian metric on D by a Finsler metric allows more general action 
functionals to be considered. The rather stringent requirement that the Lagrangian be the 
square root of a quadratic form is replaced by the more natural requirement that it be 
positive-homogeneous in u (which is necessary if the length of an oriented curve is to be 
independent of parametrization) , and that it be strongly convex (which is necessary if there 



are to exist locally minimizing paths in every direction). The problem of finding minimizing 
paths satisfying (jl.2|) is equivalent to finding geodesies of the sub-Finsler manifold (X, D, F). 

In this paper we will investigate sub-Finsler manifolds in the simplest nontrivial case: a 
three-dimensional manifold X with a rank two contact distribution D. We will work locally, 
and thus we will not generally concern ourselves with the issue of local vs. global existence 
of objects such as coordinates, vector fields, and differential forms. 

In the next two sections we will review some results of Hughen ^T] concerning sub- 
Riemannian geometry in dimension three and some results of Cartan [3 IH] concerning the 
geometry of Finsler surfaces. We will then combine these techniques to construct a complete 
set of local invariants for sub-Finsler manifolds in dimension three via Elie Cartan's method 
of equivalence. (See [S] for an exposition of this method. The reader should be aware that 
where Gardner uses left group actions, we use right group actions for greater ease of com- 
putation.) Additionally, we will derive the geodesic equations, compute the Jacobi operator 
for the second variation problem, and investigate homogeneous examples. 

2. Review of sub-Riemannian geometry of 3-manifolds 

The material in this section is taken from Keener Hughen's Ph.D. thesis ^I]. Unfortu- 
nately this thesis was never published, but some of the results are summarized in |14j . 

Let (X, D, (, )) be a sub-Riemannian structure on a 3-manifold X with a contact distribu- 
tion D. A local coframing [Tf-^rf^rf) on X is said to be ^-adapted to the sub-Riemannian 
structure ii D = {rf}^ and (, ) = {vj^Y + {''f'Y- The set of 0-adapted coframings of X forms 
a Go-structure Bq — > X, where Gq is the Lie group 

^ ^ : AG 0(2), 6GE^ cGl 
c ^ ' 
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We apply the method of equivalence to this Go-structure, and after two reductions we arrive 
at the bundle of 2-adapted coframings. This is a G2-structure !B2 -^ X, where G2 is the Lie 
group 

[A 
detA 

There is a canonical coframing {u^ , u"^ , u^ , a) (also known as an (e) -structure) on S2 whose 
structure equations are 

dto'^ = -a A uj^ + Aico'^ A uj^ + A2 cu^ A uj^ 

(2.1) du'^ = aAuo^ + A2ijo'^ Auj^ - Aiu'^ Au^ 

du^ = LU^ Au'^ 

da = Siuj'^ A uj''^ + S2Uj''^ A uj^ + K uj^ Acu^. 

Differentiating these equations shows that 

3 
dAi = -2A2a + J^ Bu ij 
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dA2 = 2Aia + ^5: 
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for some functions Bij on S2, and that 

5*1 = B12 — B21, 5*2 = Bii + i?22- 

By the general theory of (e)-structures, the functions Ai,A2,K form a complete set of 
differential invariants for the G'2-structure S2 -^ X, and hence for the sub-Riemannian 
structure (X, D, {,)). 

For later use, we observe that S2 may be viewed geometrically as a double cover of the 
unit circle bundle of the sub-Riemannian metric. If the sub-Riemannian structure (X, D, {,)) 
is orientable (i.e., if we can choose an orientation on each of the subspaces D^ which varies 
smoothly on X), then S2 consists of two disjoint connected components. In this case we 
can restrict the set of 0-adapted coframings by requiring that such a coframing be oriented, 
i.e., that the 2-form rj^ A r^^ be a positive area form on D. Doing so replaces the 0(2) 
component of the structure group by S0{2). This does not change anything essential in 
the preceding discussion, but it does lead to a G2-structure B2 which is connected and is 
naturally isomorphic to the unit circle bundle of (X, D, (, )). 

3. Review of Finsler geometry of surfaces 

The material in this section is taken from jB]. (We will, however, use the more standard 
notation for the invariants which is found in p.) 

A Finsler metric on a surface M is determined by its indicatrix bundle: this is a smooth 
hypersurface S^ C TM with the property that each fiber S^; = EflTrM is a smooth, strongly 
convex curve which surrounds the origin Oa; G T^-M. A 3-manifold S C TM satisfying this 
condition is called a Finsler structure on M. A differentiable curve 7 : [a, 6] — ^ M is called 
a Tj-curve if, for every s G [a,b], the velocity vector 7'(s) lies in S. The following result is 
taken from [3] and is due to Cartan jE]: 

Proposition . Let S C TM be a Finsler structure on an oriented surface M, with basepoint 
projection tt : E — >■ M. Then there exists a unique coframing {u^ , uj"^ , a) on S with the 
following properties: 

(1) uj'^ Auj"^ is a positive multiple of any ir-pullback of a positive 2-form on M. 

(2) The tangential lift 7' of any T,-curve satisfies {'y')*uj'^ = and {'j')*uj^ = dt. 

(3) du^ Auj'^ = duj'^ Aa = 0. 

(4) u^ Aduj'^ = J^AdiJ'. 

(5) diJ" = -a Auj"^. 

Moreover, there exist functions J, J, K on T, such that 

dco^ = —a A iiP' 
(3.1) du? = aAu^-IaAu? 

da = K uj^ Au? ^ J a Au?. 

The Finsler structure on M is Riemannian if and only if / = 0; in this case, differentiating 
()3.1|) shows that J = as well, and we recover the familiar structure equations 

duy" = —a A u'^ 

doj'^ = a Aoj^ 

da = Kix)^ AiJ' 



for an orthonormal coframing {uj^,uj'^) on M. In this case, a is the Levi-Civita connection 
form, and K is the usual Gauss curvature on the surface. For general Finsler surfaces, the 
function K (called the flag curvature) is a well-defined function only on S, not on M. 

4. The sub-Finsler equivalence problem 

Let (X, D, F) be a sub-Finsler manifold consisting of a three-dimensional manifold X, a 
rank two contact distribution D on X, and a sub-Finsler metric F on D. (Recall that D is 
contact if, for any two vector fields Vi, V2 locally spanning D, the vectors Vi, V2, and [vi, V2] 
span the tangent space of X at each point.) As in the Finsler case, the sub-Finsler metric F 
is completely determined by its indicatrix bundle 

S = {uGD I F(u) = 1}. 

S has dimension four, and each fiber S^; = S n D^ is a smooth, strongly convex curve in D^ 
which surrounds the origin 0^; G D^. A 4-manifold S C TX satisfying this condition will be 
called a sub-Finsler structure on (X, D). 

We will compute invariants for sub-Finsler structures via Cartan's method of equivalence. 
We begin by constructing a coframing on S which is nicely adapted to the sub-Finsler struc- 
ture; this procedure closely follows that used in j2] for constructing an adapted coframing 
for a Finsler structure on a surface. 

Let g be any fixed sub-Riemannian metric on (X, D), and let Si be the unit circle bundle 
for g. Then there exists a well-defined, smooth function r : Si — i> R"*" with the property that 

E = {r(u)~^u I uG Si}. 

Let p : S ^ Si be the diffeomorphism which is the inverse of the scaling map defined by r; 
i.e., p satisfies 

p(r(u)^^ u) = u 

for u G Si. 

Let TT : S ^ X, TTi : Si ^ X denote the respective basepoint projections, and let 
u G S. (We trust that using the same notation for points in S and in Si will not cause 
undue confusion.) We will say that a vector v G T^T, is monic if 7r'(u)(v) = u. Since 
7r'(u) : TuS -^ T7r(u)X is surjective with a one-dimensional kernel, the set of monic vectors 
in TuS is an affine line. A nonvanishing 1-form ^ on S will be called null if ^(v) = for all 
monic vectors v, and a 1-form a; on S will be called monic if a;(v) = 1 for all monic vectors 
V. The set of null 1-forms spans a two-dimensional subspace of T^S at each point u G S, 
and the difference of any two monic 1-forms is a null form. 

In the sub-Riemannian case, u^ is a monic form and the null 1-forms are spanned by u"^ 
and a;^. (Recall that these forms are part of the canonical coframing on Si described in 
section|21) Moreover, D is defined hy D = {u^}-^; this makes sense because according to the 
structure equations (j2.ip . u^ descends to a well-defined form on X. Since the diagram 

S^— Si 

X 

commutes, it is straightforward to verify that the null forms on S are spanned by p*{uj'^) and 
p*{uj^), that D = {p* {uj^)}-^ , and that p*{ruj^) is a monic form on S. 



A local coframing [r]^ , rj"^ , rj^ , (p) on S will be called 0-adapted if it satisfies the conditions 
that f/^ is a monic form, f/^ and f/^ are null forms, and D = {r/^}"*". For example, the coframing 
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(4.1) 7]^ = p[ruj^}, ri~ = p[uj'), ri" = p[uj"), (p = p[a) 

is 0-adapted. Any two 0-adapted coframings on S vary by a transformation of the form 
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63 



(4.2) 
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with 6163C4 7^ 0. The set of all 0-adapted coframings forms a principal fiber bundle Sq — ^ S, 
with structure group Go consisting of all matrices of the form ()4.2|) . The right action of Go 
on sections cr : S —i> Bo is given hj g ■ a = g^^cx. (This explains the inverse occurring in 

There exist canonical 1-forms r^^, 77^, 77^, on Bo with the reproducing property that for any 
local section cr : S ^ Bo, 

a*{r]')=r]\ a*(0) = 0. 
These are referred to as the semi-basic forms on Bo. A standard argument shows that there 
also exist (non-unique) 1-forms aijjSij'ji (referred to as pseudo- connection forms or, more 
succinctly, connection forms) , linearly independent from the semi-basic forms, and functions 
Tj^ on Bq (referred to as torsion functions) such that 
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These are the structure equations of the Go-structure Bq. The semi-basic forms and connec- 
tion forms together form a local coframing on Bq. 

We proceed with the method of equivalence by examining how the functions T*^ vary if 
we change from one 0-adapted coframing to another. A straightforward computation shows 
that under a transformation of the form ()4.2|) , we have 



Tin — 04^1, 
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In particular, the functions TiQ,Tf2 are relative invariants: if they vanish for any 0-adapted 
coframing, then they vanish for every 0-adapted coframing. The coframing ()4.H) has T^q = 
r~^, so we can assume that these invariants are nonzero. ()4.4|) then implies 



7-3 



that we can adapt coframings to arrange that 
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A coframing satisfying this condition will be called 1-adapted. For example, if we set 

dr = riuj^ + r2 u? + r3 u;'^ + Tq 0, 
then the coframing 

(4.5) f]" = p*{ru^ -Tquj^), ff = p*{ruj^), ff = p'ir'^iJ'), (j) = p*[a) 

is 1-adapted. Any two 1-adapted coframings on S vary by a transformation of the form 
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with b\ 7^ 0. The set of all 1-adapted coframings forms a principal fiber bundle Si C Bq, 
with structure group G\ consisting of all matrices of the form ()4.6p . When restricted to Si, 
the connection forms ai, /^s — /3i, 74 — /3i become semi-basic, thereby introducing new torsion 
terms into the structure equations of !Bi. By adding multiples of the semi-basic forms to the 
connection forms so as to absorb as much of the torsion as possible, we can arrange that the 
structure equations of Si take the form 



(4.7) 
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-T]^ A (j) 

r/^ A r/2 + Tf^ r^i a r^^ + T^^ rf A 




= d{drf) mod rf 



therefore, T|q 



0. 



We now repeat this process. Under a transformation of the form (j4.6p . we have 

-^20 ~ "1-^20 
(4.8) T/2 = ^12^/2 - biCiT^o - «2 
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In particular, T2Q is a relative invariant which transforms by a square, so its sign is fixed. 
The coframing ()4.5|) is 1-adapted, and if we set 

dro = roi u^ + ro2 uj'^ + roa uj^ + roo 0, 

it has T2Q = . The condition that each fiber of S be a strongly convex curve enclosing 

the origin is exactly the condition that this quantity be positive (see Lemma I?. 31 for a proof), 
so we can assume that T2Q > 0. (j4.8|) then implies that we can adapt coframings to arrange 
that 
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A cofraining satisfying this condition will be called 2- adapted. Any two 2- adapted coframings 
on S vary by a transformation of the form 
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with e = ±1. The set of all 2-adapted coframings forms a principal fiber bundle S2 C Bi, 
with structure group G2 consisting of all matrices of the form ()4.9|) . When restricted to S2, 
the connection forms Pi, 71 + 0:2, 72 + 2/32 become semi-basic. By adding multiples of the 
semi-basic forms to the connection forms so as to absorb as much of the torsion as possible, 
we can arrange that the structure equations of B2 take the form 
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Moreover, we have 



?7^ A 
-r/i A </) + Tfa r]^ Ar]"^ + T^q rf A (f) 
T]^ At]"^ + Ti^ 7f Ar]^ - T^2 T]^ At]^ + T^q t]^ A 
T°2 r^i A r/2 + T°o r/^ A + T°o r/^ A 
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Under a transformation of the form ()4.9|) . we have 
(4.11) f22 = T22 + £(a2r|o + &2) 



^23 — ^^2.3 + 2&2^c 
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A coframing satisfying this condition will be called 3- adapted. Any two 3-adapted coframings 
on S vary by a transformation of the form 
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The set of all 3-adapted cofrainings forms a principal fiber bundle IBs C B2, with structure 
group 6*3 consisting of all matrices of the form (J4.12J) . When restricted to S3, the connec- 
tion forms 02, (^2 become semi-basic. By adding multiples of the semi-basic forms to the 
connection forms so as to absorb as much of the torsion as possible, we can arrange that the 
structure equations of S3 take the form 
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7]^ A + T^r^ r^i A ?7^ + T2^3 7]'^ Ar]^ + T^^ rp Acf) 

-77^ A + T^Q rf A(j) + Tfg 77^ A 77^ + T^^ rf At]'^ + T^q t]^ A > 

7]^ Arf + T^^iff A 

TO2 r/ 1 A 77 2 - T3I0 r/ 1 A + T^o r/ 2 A 

(The coefficients Tj'g, T20 in (J4.13J) are slightly modified from those in (J4.10J) .) Moreover, we 
have 

= d{drf) mod 

= -(T/3 + T|3 + T2W2)^^Ar7^Ar/^ 



therefore, we can set 
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for some function A2 on S3. (The reason for this choice of notation will shortly become 
apparent.) 

Under a transformation of the form (|4.12p , we have 

(4.14) f^,=T,\ + ecs 
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so we can adapt coframings to arrange that 

rpl 

-^23 

for some function Ai. A coframing satisfying this condition will be called 4-O'dapted. Any 
two 4-adapted coframings on E vary by a transformation of the form 
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The set of all 4-adapted coframings forms a principal fiber bundle S4 C B3, with structure 
group G4 = Z/2Z. B4 is thus a double cover of S, and the 1-forms {1]^, rf, if, (p) form a 



10 



canonical coframing on S4. When restricted to B4, the last remaining connection form 73 

becomes semi-basic, and the structure equations of B4 take the form 

(4.16) 

r^^ A - (A2 + irloTOa) r/i A 7^3 + Ai 7^2 A r/3 + T3I0 7^3 A 



drf 
drf 
dcf) 



-T]^ A(P + T^o r/2 A + Ai r^i A r/3 + (A2 - ^T^f^T^^) rf Arj^ + T^^ rj^ A 



ri^ Arf + T^qT]^ A 

T^2 ^^ A r/2 + T°3 7]^ Ar]^ + T°3 rf Ar]^ - T^q t]^ A (f) + T^q rf A (j) + T^q r]^ A i 

Finally, we differentiate equations ()4.16p in order to find any remaining relations among the 
torsion functions. Setting 

and computing d{drf) = yields 

(4.17) ^2^1 = ^10 + 7^10^30 

rp2 (rpl I rp2 rpQ \ 

-^ 20,2 ~ l-^ 30 "T -^ 20-^ 20 J • 

Then computing d{drf) = mod rf yields 



rpQ 

-^20 



27^2 
^-^30- 



Further differentiation yields only differential equations for the torsion functions and no new 
functional relations. 



If we rename the T*^ as follows: 



rpl 
-'20 

-^30 

rp2 

-^30 

rpO 

-^12 

rpO 
-'30 
rpQ 
-'23 
rpO 
-'13 



/ 

Jl 
J2 

K 

So 
Si 
-S2, 



drf 
drf 
drf 



then the structure equations on S4 become 
(4.18) 

rf A(j) + AiT]'^ Ar]^ + (^2 + \IK) r]^ A r]^ + Ji r]^ A (j) 

-1]^ A + (A2 - ^IK) rf A rj'^ - Airf A 7]^ + J2rf A (p + I rf A (p 

rf Arf + I rf A(f) 

Sorf A(j) + Sirf Arj'^ + S2rf Arf - Jirf A(l)-2J2rf A(l) + Krf A rf 

(Compare with the sub-Riemannian structure equations ()2.H1 .) 

Our first result is that I is the fundamental invariant that determines whether or not a 
sub-Finsler structure is sub-Riemannian: 

Theorem 4.1. The sub-Finsler structure S is the unit circle bundle for a sub-Riemannian 
metric if and only if I = 0. 
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Proof. One direction is trivial: if S = Si for some sub-Riemannian metric, then the canonical 
coframing which we have constructed on S is simply 

{r]\ if, if, 0) = (a;\ u^, u^, a), 

and so the structure equations (|4.18p must reduce to (|2.1|) : therefore, 7 = 0. 
Now suppose that 1 = 0. Then 

= d{dri^) = ( Ji r/2 - J2 T]^) A r/3 A </.; 
therefore, Ji = J2 = 0. Now computing d{dr]^) = mod rj^ shows that 

dAi = {So - 2A2) mod 7]\ rf , rp , 

and computing d{diif) = mod r/^ shows that 

dAi = {-So - 2A2) (p mod r/^ r/^ i]^. 

Therefore, 5*0 = 0, and the structure equations ()4.18|) have the form ()2.1|) . This implies that 
S is the unit circle bundle for a sub-Riemannian metric, as desired. D 

5. The geodesic equations 

In this section we consider the problem of finding geodesies of the sub-Finsler structure. 
Recall that the sub-Finsler length of a horizontal curve 7 : [a, 6] — i> X is given by 

(5.1) £(7)= f'F{i{t))dt. 



Finding critical points of this functional amounts to solving a constrained variational prob- 
lem. However, care must be taken when computing variations among horizontal curves on a 
non-integrable rank s distribution D. Given a horizontal curve 7, one would like to consider 
"D-variational vector fields on 7 that vanish at the endpoints," but in general the existence 
of such vector fields is far from guaranteed. In fact, this can fail spectacularly: for example, 
when D is an Engel system on a 4-manifold M, M is foliated by horizontal curves that have 
no such variations f3^. 

If such a vector field exists along 7, then 7 is said to be regular, and the methods outlined 
in Pj suffice to find the first variation. A horizontal curve for which this fails is called non- 
regular (or abnormal). In |10' Lucas Hsu established the following criterion for a curve to 
be non-regular: 

Theorem 5.1. (Hsu, ^Dlj LetX C T*X be the annihilator of the ranks distribution D C TX, 
and let \I' be the pullback of the canonical symplectic 2-form on T*X to X. A horizontal curve 
7 : [a, 6] — >■ X is non-regular if and only if it has a lifting 7 : [a, 6] — >■ X that does not intersect 
the zero section and satisfies 7'(t) _i\E' = for all t E [a, b]. 

In the present case, D is a contact system on a 3-manifold with X = span {rj^}, and it is 
easy to see that in this case all horizontal curves must be regular. In what follows we will 
therefore use the variational methods described in [H]; our argument closely follows that of 

Choose an orientation of D, and consider the set of coframes in S4 that preserve this 
orientation; for simplicity we will continue to use the notation B4 for this set. Every hori- 
zontal curve 7 : [a,b] —>■ X lifts to an integral curve 7 : [a, b] —>■ B4 of the differential system 
X = {rj^,rj^} with rj^{'y'{t)) 7^ 0. This lift corresponds to choosing a 4-adapted coframing 

along the horizontal curve so that the vector ei dual to rj^ points in the direction of the 
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velocity vector of the curve. The sub-Finsler length of 7 is then equal to the integral of the 
monic one- form 7]^ along the lifted curve 7. The problem of finding critical curves of the 
sub-Finsler length functional among horizontal curves is thus equivalent to finding critical 
curves of 



(5.2) £(7) = / v' 

among integral curves 7 of X = {rj'^, rj^} on B4. 

Proposition 5.2. The critical curves of C among integral curves ofT on S4 are precisely 
the projections of integral curves, with transversality condition rj^ 7^ 0, of the differential 
system J = {rf^rj^^ip — Xri^,dX — Ct]^} on V = B4 x M, where A is the coordinate on the 
fiber R and C = A^J + A Ji + A2 + \IK. 

Proof. Following the algorithm in jH], we define a submanifold 2, C T*!B4 as follows: for each 

X G S4, let Zx = 77^ (x) + span{Xr} and let 

(5.3) Z = U Z,. 

xe'B4 
Let C be the pullback to Z of the canonical 1-form on T*^)^. By the "self-reproducing" 
property of C, we may write 

(5.4) C = ^^ + ^2V'^ + hv^ 

(where we have suppressed the obvious puUbacks in our notation). According to the general 
theory described in jH] , the critical points of the functional 



(5.5) C{l)= IC 

among unconstrained curves 7 on Z project to critical curves of C among integral curves 7 
of X on B4; moreover, a curve 7 on Z is a critical curve of C if and only if 7'(t) _i(iCl7(t) = 0. 
A straightforward computation shows that 

rfC = A2 A r/^ - (1 + A2/)0 Av^- (A3/ + ^1 + A2 ^2)0 A 7]^ + X^T]^ A r/2 

(5.6) + {A2 + \IK - X2A{)Tf' Ari^ + iA^ - \X2IK + X2A2W A v"" 

+ dX2 Ar{^ + rfAs A rf". 

By contracting dC, with the vector fields dual to the coframing {r^^, rf.rf, 0, dX2, dX^} on Z, 
we find that subject to the condition 7*77^ 7^ 0, the requirement that 7'_i(i(^ = is equivalent 
to the condition that 7 is an integral curve of the system 

J = {ri\ ri\ - X^ri\ dX^ - {Xp + A3 Ji + A2 + \lK)r^^} 

on the submanifold V C 2. defined by A2 = 0. (Henceforth we will omit the subscript on 
A3.) Curves satisfying this requirement project to critical curves of the functional £ among 
integral curves of X on B4, and thus to local minimizers of the sub-Finsler length functional 
C on X. Since every horizontal curve on X is regular, every local minimizer arises in this 
way. 

D 
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We will call a unit speed horizontal curve 7 : [a, 6] — > M a sub-Finsler geodesic if it has a 
lift to an integral curve of J7 on "y. When 7 has unit speed, it lifts to an integral curve of JT 
if and only if it satisfies the geodesic equations 

(5.7) 7]^ = ds, 1]'^ = 0, 7]^ = 0, (p = \ ds, d\ = C ds. 

6. The Jacobi operator and the second variation 

This argument is similar to that given in JT] for the sub-Riemannian case; we will describe 
it in some detail since |TTj is unpubhshed. 

Since the geodesic equations are defined on the bundle V = B4 x M, we will work on ^ and 
use the coframing {77^, 77^, r]'^,-?]^, 77^}, where 

?7^ = — Xr]^ 



V 
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d\-{\^I + \Ji + A2 + yK)r]\ 



The structure equations ()4.18|) imply that this coframing has structure equations 

dri'' = r]^ A T]^ - Xr]^ A T]"^ + AiT]'^ A T]^ + JiT]^ A T]^ - (A Ji + A2 + \IK) r]^ A rj^ 
drf = -7]^ Ar]^ + J2rf At]^ + 17]^ At]^ -Xlr]^ Arf 

+ {A2 - \IK) r]^Ar]'' + {Ai - A J2) r]' A if 

(6.1) drf = 7]^ Arf + Irf At]^ -Xlr]^ Arf 

d7]^ = r^i A r/^ - Ji r^i A r/^ + {Si - XAi) rf Arf + {Sq - AJi) jf A r]^ 

- (A + 2J2) rf Ati^ + (A^ + 2XJ2 + K)r]'A if 

+ (JiA^ + {A2 - So + \IK)X - S2W A V^ 
drf = -(A^J + A Ji + A2 + \IK) [if Aj]'^ - Xt]^ Aif + Anf At]^ + Jiff A r]^ 

- (A Ji + A2 + \IK) 7]' A f]^] 

+ 7]^ Ad{X^I + XJi + A2 + \lK). 

As in the previous section, every horizontal curve 7 has a canonical lift to an integral 
curve 7 of the system X = {rj^, rf} on B4. This in turn has a canonical lift to an integral 
curve 7 of the system X = {77^, rf^ r^^} on y. The length of 7 is equal to the integral of rf 
along the lifted curve 7, and 7 is a geodesic if and only if 7 is an integral curve of the system 
J = {rf^ jf^ 774 7^5j Qj^ \^ 

Suppose that 7 : [0,^] -^ X is a horizontal curve joining points p and q in X. If 7^ is a 
fixed-endpoint variation of 7 through horizontal curves, then 7^ lifts to a variation 7^ of 7 
through integral curves of X; this variation does not necessarily fix endpoints, but it satisfies 
the condition 

vr 07^(0) =79, 7ro74(£) = g, 

where vr : y ^ X is the usual base point projection. A variation 7^ satisfying these conditions 
will be called an admissible variation of 7, and its variational vector field -^ at t = will 
be called an infinitesimal admissible variation along 7. 

Now suppose that 7 is a geodesic. Let 7t_„ be 2-parameter fixed-endpoint variation of 7 
through horizontal curves, and let '^t,u be its lift to V, with infinitesimal admissible variations 

d ^ d ^ 

^(■5) = ^^*."('^) \t=u=0, W{s) = -^1t,u{s) \t=u=0 ■ 
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Let (ci, . . . , 65) be the framing dual to the coframing [t]^, . . . , rj^) on y, and write 

5 5 

V{s) =J2Vi{s)e,{s), W{s) = Y,W.{s)ei{s). 

i=l 4=1 

The Hessian £**(V^, W) of the length functional is, by definition, 

C^^{V,W) = Q^^ilt,u) \t=U=0 ■ 

Proposition 6.1. Let 7 : [0, £] — ^ ^ 6e a lifted geodesic with 2-parameter admissible variation 
7t^„. The Hessian of the length functional evaluated at the infinitesimal admissible variations 
V = ^ ViCi and W = J2 ^i^i ^^ 

C,,{V,W)= [ W3JiV3)ds, 
Jo 

where J is a self-adjoint, fourth-order differential operator on the space of smooth functions 

on [0, i] given by 

(6.2) J{u) = -u + — (Pii) + Qu, 

as 

for certain functions P, Q along 7. 

Here the dots over u represent derivatives with respect to s, and the precise definitions of 
P and Q will appear in the proof. 

Proof. For any smooth function / on V, we will write 

5 

df = J2l^'^'- 

i=l 

Differentiating (j6.ip yields relations among the derivatives of the invariants of the sub-Finsler 
structure, and these relations must be taken into account in the computations that follow. 
As in illj, the Hessian £**(y, W) is equal to the integral 

W_\d{V^dC), 

7 
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where ( = rj^ + Xr]^. A long but straightforward computation shows that along 7, the 
integrand W ^d{V ^dQ is equivalent modulo J' to 

W2 (V^ + A1V3- (A' + 2 J2A + Ai + K)V2 

+ ( - JiA^ + {AJ - \IK -A2 + Sq)\ + S2)V^ + JiVi - V^ 

+ W^3 ( - V^5 - ^iV^2 + (/A + Ji)\/4 

+ ( - JiA' + (All - \IK - A2 + So) A + ^2 - Ai,i)V^2 

+ ((So/ + 1,3 + So,4)A2 + ((Ai + K)J2 + JlSo - Si - So,l + S2,4)A 
(6.3) + {\PK^ + i(//^,3 + /s/^) + A? + ^2 + A2IK + J1S2 + A2,3))^3 

+ (/4A2 + (JJi + Ji,4)A + \{PK + ir/4 - Mi,4) 

- Ai/2 - A2I - 2SqI + J2 - Ai + A2,4)V^4 

+ /Al^s) 

^( - \/2 - (/A + Ji)V^ + Ji^s + (-/'A^ - {IJ^ + J2)A + ^1)^3 - V^ 



+ Wa 



(This computation takes into account the fact that along 7, Vi^i = V^.) 

Now define T(t,u,s) = 7t,u(s). Since each curve ^t,u is an integral curve of T, we have 

Vi{t, u, s) dt + Wiit.u, s) du + Yi{t, u, s) ds 
V2{t, u, s) dt + W2(t, u, s) du 
T* v^ = V3{t,u,s)dt + W3{t,u,s)du 
¥4(1, u, s) dt + W^it, u, s) du 
V^it, u, s) dt + W^it, u, s) du + Y^it, u, s) ds 

for some functions Vi,Wi,Yi satisfying (with a minor abuse of notation) Vi{0,0,s) = Vi{s) 
Wi{0, 0, s) = Wi{s). Since 7 is the lift of a geodesic, we also have 



'ri'' 




v' 




v' 


= 


v' 




v' 





YiiO,0,s) 



1, 



n(0,0,s)=0, 



dYi 



(0,0, s) 



dYi 



(0,0,s) = 0. 



dt ^"' "' ' du 

When the structure equations (j6.1|) are pulled back by T and then restricted to 7, they imply 
that 

Vi = -XV2 - (JiA + A2 + \IK)V3 

(6.4) V2 = -I\V2 + {-J2\ + A^)V3-V4 

V3 = V2- IWs 

n = (A' + 2J2A + K)V2 + (JiA^ + {A2 + \IK - So)A - ^2)^3 - J1F4 + ^5- 

The third equation in ()6.4j) implies that V2 = V3 + /AV3. The first equation in ()6.4|1 can then 
be written as 

1(1-, + AV3) = 0. 
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Since V is an admissible infinitesimal variation, V^i, V2, and V3 must vanish at the endpoints 
of 7, and it follows that Vi = — AV3. Equations (j6.4j) can now be used to express Vi, V2, V4, 
and V5 in terms of V3 and its derivatives, as follows: 

Vi = -V3 - 2IXV3 + [-(2/2 + /4)A2 - (2/Ji + 2J2)A + (Ai - IA2 - IPK)]V3 
(6.5) V, = -V3-{2IX + Ji)V3 

+ [-(4/2 + 3/4 + 1)A2 - (8/Ji + 6J2)A + (Ai - 3/A2 - |/2/C - KM 
+ [-(4/=^ + 6//4 + /)A3 

- (12/Vi + 5/4J1 + /41 + 2/Ji,4 + 8/J2 + 2J2,4 + ^i)A2 

+ (-3//S:/4 - 3A2/4 + 2Ai,4 + /^Ai,4 - 2/^2,4 + 4Ai/ + 2/=^Ai 

+ 3^2 - 21^ K - MJl - 8J1J2 - pK + 3So + APSo)X 
+ (Ai,i - /A2,i - |/2/^,i + AiJ, - AIA2J1 - 3A2J2 - IPJiK 

- 2IJ2K + S2)]V3. 

When these expressions are substituted into the integrand ()6.3p . it takes the form 

Ws{V3 + f^{PV3) + QVs), 
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where 



P=i2I^+4l4+l)X'^+8(JiI+J2)X+[2KI^+4A2l+K-3Aij 
Q=[4/4+2(8/4+1)/2+{5/24+/4)]A4 



4" 
+ [20Ji/3 + (14J2 + 10Ji,4)/2 + (40Ji/4+6/41+4Ji+8J2,4)/+(12J2/4+6Ji,4/4+<^2+Ji.4)]A3 

+ [(3E'-8Ai)/4-(ioA2+165o+4yli,4)/^ + (16E'/4-10Ai/4-19Ai+8A2,4+40Jf+ii')/2 

+ {5/12/ 4-4yll,4/4-165o/4+42JiJ2 + 18JlJl,4+2Ji,4l-12A2-^^l, 4-650)/ 
+ {5/i'/24 + 18j2/4-13Ai/4+6A2,4/4+2/^/4+6Jl/41+/411+/3-2Ai+A2,4 
+3j2+8j|+8J2Jl,4 + 10JlJ2,4+2J2,41+5'o,4)]A2 

+ [{12/s'Ji-20AiJi-2Ai,i+ft:,i)/3 

+{4/i'Ji,4-9JiAi,4-2A2,i-Ai,4i-14JiA2-12J2Ai+9J2/i'-36JiSo-45o,i)/2 

+ {20/^ Jl/ 4+4/^,1/ 4+4/^/ 41+7A2Jl,4+4/i'J2,4 + 16JlA2,4-5J2Ai,4-3Al,i+2A2,41 

+ |/^,i-34AiJi-4A2J2+4/^Ji+24j3-20J25o+2S2)/ 

+ {14JiA2/4+4A2,l/4+8/^J2/,4+5A2/41+7A2J2,4"8JlAi,4+8J2A2,4-3A2,l-2Al,41 

-4S'o,i+S2,4-18AiJ2-13A2Ji+24j2j2+|/<'J2-8Ji5o-S'i)]A 

+{^K'2Ia-Ia2Ai^4+3KA2,4+3JiK^i + ^K^ii~7A1+K'^-YAiK+9KJ'1~10A2So)P 
+(6A2/^/4+5A2A2,4+4JiA2,i-3/^Ai,4+A2,ii+3J2/i',i + |/^,3-llAiA2+12A2Ji' 

-3A2/i'+llJlJ2/i'-|/^S'o)/ 

+{3A2/4+i/^/3-5A2Ai,4+4J2A2,i+A2,3-Ai, 11-52, 1+2A2-8A2+12A2J1J2 
-6A25o+2j|/^+Ji52)]. 

Thus the proposition is proved. D 

We saw in the proof of this proposition that any infinitesimal admissible variation V = 
T^iViCi satisfies 

and that Vi, V2, V3 vanish at the endpoints 7(0), 'j{i) of 7. In particular, V3 vanishes to first 
order at and i. Let C^[0,£] denote the space of smooth functions on [0,i] that vanish to 
first order at the endpoints, and note that the Jacobi operator J is formally self-adjoint on 

Define a quadratic form Q{u) by 

Q{u) = C^^{u,u). 

Recall that the index of Q is the dimension of the largest subspace of C^[0,£] on which Q 
is negative definite. Because J is self-adjoint on C^[0,£], its eigenvalues form a countable 
subset of the real numbers with -|-oo as the only possible cluster point. It follows that J has 
only finitely many negative eigenvalues, and that therefore the index of Q is finite. 

Definition 6.2. A point c G (0, i) is a conjugate point of J with multiplicity m if the space 
of solutions of the system 

(6.6) J{u) = 0, m(0) = m(0) = 
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which vanish to first order at c has dimension m > 0. The point 7(c) along a geodesic 7 is 
a conjugate point of '-/ if c is a conjugate point of the corresponding Jacobi operator J along 

7- 

Note that, since J is a fourth-order operator, the multiphcity of any conjugate point of J 
is either one or two. 

Theorem 6.3. The index of Q is equal to the number of conjugate points of J, counted with 
multiplicity. 

Proof. Suppose that the index of Q is n, and let 

Ai < A2 < • ■ ■ < A„ < 

be the negative eigenvalues of J. For each s G (0,£], let 

Ai{s) < A2{s) < ■ ■ ■ 

denote the eigenvalues of the operator J on the space [0, s]. (Note that Aj(£) = Aj for 
i = l,...,n.) It follows from general theory (see, e.g., ^2]) that each Aj(s) is a strictly 
decreasing, continuous function on (0, i], with lim5^o+ Aj(s) = +cxo. Therefore, each function 
Aj(s), i = 1, . . . ,n has exactly one root Cj. 
These roots 

< ci < 02 < ■ ■■ < On < i 
are precisely the conjugate points of J between and i. To see this, note that by definition, 
the condition Aj(q) = implies that there exists a function Ui G C^[0, q] satisfying J{ui) = 
0. By extending Ui to a solution of J{u) = on C^[0,i], we get a solution to ()6.6|) which 
vanishes to first order at q. Conversely, if c G (0,£) is a conjugate point, then J has a zero 
eigenvalue on C^[0, q]; therefore, Aj(c) = for some positive integer j. Since A^ is a strictly 
decreasing function, it follows that Aj{i) is equal to one of the negative eigenvalues Xj of J, 
and therefore, c = Cj. D 

Corollary 6.4. A geodesic no longer minimizes length beyond its first conjugate point. 

Proof. Suppose that 7(c) is the first conjugate point of the geodesic 7. Theorem 16 . 31 implies 
that for every i > c, the index of Q on the space C^[0,£] is positive, and so there exists a 
function u G C^[0, i] for which Q{u) < 0. Setting V3 = u and defining Vi, V2, V4, V5 according 
to equations (j6.5|) defines a direction V = J2i=i ^i^i along which the length functional C (and 
hence C) decreases. D 

7. Symmetries and homogeneous examples 

In this section we examine the symmetries of sub-Finsler structures and describe homo- 
geneous examples. 

Definition 7.1. Let H be a sub-Finsler structure on {X,D). A symmetry 0/ S is a diffeo- 
morphism $ : X ^ X which satisfies $'(S) = S. (Note that this implies that $'(D) = D 
as well.) A symmetry of the 'L/2'L- structure B4 is a diffeomorphism \1/ : B4 — > S4 with the 
property that 

(^^*r]\ ^*r/2, ^V, ^*0) = (r^i, r^^, r]^, 0). 
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A standard argument shows that the map 

gives a one-to-one correspondence between orientation-preserving symmetries of S and sym- 
metries of B4. By a theorem of Kobayashi ^^l, it follows that the group of symmetries of S 
can be given the structure of a Lie group of dimension at most four. 

There are two possible definitions of homogeneity for a sub-Finsler structure: we could say 
that S C TX is homogeneous if its group of symmetries acts transitively on X, or we could 
require the more restrictive condition that this group act transitively on S. Both notions 
are interesting, and we will consider each of them in the remainder of this section. 

7.1. Symmetry groups of dimension four. First we consider the case where the group 
of symmetries of S is four- dimensional and acts transitively on B4. Since any symmetry 
must preserve the canonical coframing [i]^, 77^, 77^, 0) on B4, it follows that all the torsion 
functions must be constants. Conversely, if all the torsion coefficients are constants, then the 
structure equations of S4 define a local Lie group structure on S4 for which the canonical 
coframing is left-invariant; this Lie group then acts transitively on S4 in the obvious way. 

So, suppose that all the torsion functions in the structure equations ()4.18|) are constant. 
Then 

= d{dr]^) = [-(/Ji + J2) v^ + (Ji - 2/J2) 7]^] Ar]''A(P; 
therefore, Ji = J2 = 0. Next we have 

= d{dr]^) = -S2 if- A Iff A if mod 
= d{drf) = Siif Arf A if mod 0; 
therefore, ^i = 5*2 = 0. Then 

= d{d<i)) = {So - IK) if A if A <p mod rf; 
therefore, Sq = IK. Now 

= d{di]') = [(2Ai + IA2 + \PK) 1]' + {-21 Ai + 2A2 - IK) 1]^] A r/3 A 0; 

therefore, 

PK _ IK{P-2) 

'~~7^T2' '~ 2(/2 + 2) ■ 

Finally, we have 

Q = d{dif) = j:^-^ri^ Ari^'Act); 

therefore, IK = 0. If / = 0, then S is sub-Riemannian; the homogeneous sub-Riemannian 
structures are classified in [TTj. So suppose that / 7^ 0. Then we have K = Q, and the 
structure equations ()4.18j) reduce to 

drf = if /\(j) 

(7.1) dif = -1]^ A (j) + I if A (j) 

drj-^ = if A rf + I if A (I) 

d4> = 0. 

Differentiating ()7.1|) yields no additional restrictions; thus there exists (at least locally) a 
1-parameter family of homogeneous sub-Finsler structures which are not sub-Riemannian. 
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In fact, these structure equations can be integrated explicitly. First, since d(f) = 0, there 
exists a function on S such that 

(I) = d9. 

Next, note that the system S = {t]^, t]"^} is Frobenius (i.e., dS = mod S); therefore, there 
exist functions x, y on E, independent from 9, such that 

S = {dx, dy}, 

and functions aij, i,j = 1,2, such that 

1]^ = ail dx + ai2 dy 

rf = a2i dx + 022 dy. 
Now the first two equations of ()7.H) imply that the aij are functions of x, y, 6 alone, and that 
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are each sc 






df 
86 


= -9 












dg 

06 


= f- 


-I9 







(7.2) 



for functions /(x, y, 6), g{x, y, 6). The solution of these differential equations depends on the 
value of /. 

7.1.1. Case 1: P > 4. The general solution of ()7.2|) in this case is 

f = c,ix,y)e'^' + C2ix,y)e^^' 

g = -ci{x, yyie''^^ - C2{x, y)r2e''^\ 
where 

n,r2 = i(-/±V/2-4). 

By modifying x and y if necessary, we can assume that 

an = ci(x, |/)e'^i^ 012 = C2(x, y)e''^'^ 

^21 = -ci{x,y)rie'^'^ 022 = -C2(x,y)r2e'^2^ 

Then the first two equations of ()7.1|1 imply that ci is a function of x alone and C2 is a function 
of y alone. 

Now the third equation of ()7.1|1 implies that 



drf = Ci{x)c2{y)e ^^\JP — Adx A dy mod t]^. 
By the Pfaff theorem, there exists a function z on S, independent from x, y, and 6, such that 



1]^ = Ci{x)c2{y)e ^^\/P - A[dz+ \{xdy -ydx)) . 
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Finally, the third equation of ()7.H) now implies that ci and C2 are in fact constants. Without 
loss of generality, we may assume that Ci = C2 = 1, and that our coframing has the form 

r^i = e^i^ dx + e'^^ dy 

rf = —Tie^^^ dx — r2e^^^ dy 

e-^Oy^r^l [dz + \{xdy- ydx)) 

de. 



v' 



7.1.2. Case 2: P < A. The general solution of ()7.2|) in this case is 

/ = e~^^'^[ci(x, y) cos{r6) + C2(x, y) sin(r6')] 
g = |e~™'^[ci(x, y)(/cos(r6') +rsin(r6')) 

+ C2(x, y){I sm{r9) — r cos(r6'))], 

where 



r = 171^^72. 



A similar argument to that given above shows that that we can take our coframing to be 

1]^ = e^^^''^[cos{r9) dx + sin(r0) dy] 



rf = ie [(-^ cos(r0) + r sin(r^)) dx + (/ sin(r^) — r cos(r6')) dy] 
if = — re~™ [dz + ^{xdy — ydx)) 
(I) = d9. 

7.1.3. Case 3: 1 = 2. The general solution of ()7.2|) in this case is 

/ = e-*[-ci(x, y) t + C2{x, y){l + t)] 
g = e"*[ci(x,?/)(l-t) + C2ix,y)t]. 
A similar argument to that given above shows that that we can take our coframing to be 

r^i = e~%{l + 6) dx - e dy] 

r]^ = e-%e dx + {1 - 0) dy] 

if = e^^^ {dz + ^{xdy — y dx)) 

4> = de. 

7.1.4. Case 4-' I = —2. The general solution of ()7.2|) in this case is 

f = e^[-ci{x,y) t + C2{x,y){l - t)] 

g = e^[ci{x,y){l + t) + C2{x,y)t]. 

A similar argument to that given above shows that that we can take our coframing to be 

1]^ = e^[{l-e)dx-edy] 

r]^ = e'^[edx + {l + e)dy] 

rf = e^^ [dz + ^{xdy — ydx)) 

(l) = de. 

Note that none of these four coframings have coordinate expressions which are periodic in 
the 9 variable. Consequently, in all four cases the indicatrix for the sub-Finsler metric in each 
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tangent space fails to be a closed curve. (In fact, these indicatrices are not even connected.) 
Therefore, these sub-Finsler structures exist "micro-locally" - that is, in a neighborhood of 
each point in S - but not locally. In other words, there is no open set 11 C X for which the 
sub-Finsler metric is defined on all of Tli. This is consistent with a theorem of Rund (see 
PP) which implies that for any Minkowski norm on a plane D^, the average value of / over 
the indicatrix must be zero; therefore, if / is any nonzero constant, the indicatrix cannot 
possibly be a closed, strongly convex curve in D^. 

7.2. Symmetry groups of dimension three. Now we consider the more inclusive case 
where the group G of symmetries of S is three-dimensional and acts transitively on X. Since 
S is invariant under the action of G, it is completely determined by the fiber E^^. at any point 
x G X. Conversely, if we fix a point a; G X and choose any smooth curve F C Dr^ which is 
strongly convex and encloses the origin, then there exists a unique sub-Finsler structure S 
on (X, D) which is invariant under the action of G and satisfies S^; = F. 

Thus we conclude that the sub-Finsler structures of this type are generated by choosing 
a three-dimensional Lie group G (or a quotient thereof by a discrete subgroup), a 2-plane 
D C TeG which is not a Lie subalgebra (so that it is bracket-generating), and a smooth 
curve F in D which is strongly convex and surrounds the origin. 

Example 7.2. Let J-C be the Heisenberg group, defined by 



1 y 

1 


z + lxy 

X 





1 



and let the contact structure on % be the rank two distribution 
(7.3) D = {dz + -{xdy-ydx)y. 

The existence of this global coordinate system makes it easy to describe sub-Finsler 
geodesies within the Heisenberg group. Moreover, this example is prototypical: by a theorem 
of Pfaff (see |4J, any contact 3- manifold has local coordinates (x, y, z) for which the contact 
system is given by the symmetric normal form ()7.3p above. 

We can define a homogeneous, flat sub-Riemannian metric on (!K, D) by declaring the 

vectors 

d y d ,, d X d 

Vi = h - — , V2 = 

dx 2 dz dy 2 dz 

to be orthonormal. Let Si be the unit circle bundle for this sub-Riemannian structure on 

!K, and define a coordinate 6 on Ei by the condition that, for u G Si, 

u= (cos^)Vi + (sin^)y2. 
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It is straightforward to check that V^i,V2 are left-invariant, horizontal vector fields on 
(!K, D), and that therefore any scaling function r{6) which depends on 6 alone defines a ho- 
mogeneous sub-Finsler structure on !K. It is also straightforward to check that the coframing 

f]^ = p* {{r cos ~ r' sin 6) dx — (r sin 6 + r' cos 6) dy) 
ff = p* [ \/r{r + r")[(sin6') dx + (cos6') dy] 



7^3 = p* fr^/'^y/r + r"[dz + l{xdy -ydx)]j 



P*\^^de 



on the sub-Finsler structure S defined by r{9) is 4- adapted. The invariants for this coframing 
are 

1 {rr'" + 3rV" + Arr') 
- ~2 ^(r + r"f/^ ' 
K = A, = A2 = Ji = J2 = So = Si = S2 = 0. 
In this case, the geodesic equations ()5.7j) can be written as 

cos^ 
r{9{s)) 

sme{s) 
dy = — T^rrras 
r{e{s)) 

X (s) sin e{s)+y{s) COS e{s) 

^^•') ^^ = Mm) ^^ 



dO^^^^^Xis)ds 



dX = IX\s)ds. 

Since the Lie algebra of J{ is solvable, it is no surprise to find that these equations can be 
solved by quadrature. First we prove a straightforward but useful lemma. 

Lemma 7.3. The expression r{r + r") is positive and bounded away from zero. 

Proof. Let x be any point in "K, and let S^ C D^ be the indicatrix in D^ corresponding to a 
sub-Finsler metric F. The strong convexity of F is equivalent to the condition that, for any 
parametrization {u{t),v{t)) of S^., 

(7.5) «V-nV; ^ Q 

u'v — uv' 
everywhere on Tj^ (see P^ for a proof). In particular, for the parametrization 

u{e) = R{e) cos e, v{e) = R{e) sin e, 

this is equivalent to saying that 

(7 6) r + r" _ R^ + 2{R'f - RR" ^ ^ 
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everywhere on S^., and hence r{r + r") > as weU. (Recall that r is the reciprocal of the 
radial position function R defining S.) Since S^ is compact, this quantity is bounded away 
from zero. By the homogeneity of "K, this bound is the same at every point x G CK. D 

This lemma and the geodesic equations ()7.4j) imply the following result. 

Theorem 7.4. For any homogeneous suh-Finsler metric F on the Heisenberg group !K, the 
sub-Finsler geodesies are straight lines parallel to the xy-plane or liftings of simple closed 
curves in the xy-plane. In the latter case, the simple closed curves are the curves of minimal 
Finsler arc length enclosing a given Euclidean area in the plane. 

Proof. A curve 7 : [a, 6] ^ J{ is a geodesic if and only if it is an integral curve of the system 
()7.4|1 . The equations for d9 and d\ in ()7.4|1 allow us to write 



dX ^ JT + r" ,^ 
dO 



A \ r 



2 r{r + r") 

1 d{r{r + r")) dr 



d9 



2 r{r + r") 
and so 

A 



cAo 



r^r{r + r") 

where Aq, c are constants, with c = r(0)A/r(0)(r(0) +r"(0)) > 0. 

Integral curves of (j7.4j) corresponding to Aq = are straight lines parallel to the a:?/-plane. 
If 7 is an integral curve corresponding to some Aq 7^ 0, then we have 

dQ = — -ds. 

ryr + r") 

By the preceding lemma, the quantity in the denominator is positive and bounded away from 
zero; thus 9 varies monotonically with s, without bound. We may therefore reparametrize 
the equations for dx, dy and dz in terms of 6. If (xq, yo) = {.x{6q), y{9o)) is any point on the 
projection of 7 to the xy-plane, then for any other value 6, we have 

1 [^ 
(7.7) x(^) - xo = ^ / (r + r") costdt 



cAo J 00 



1 /-^ 

y{0)-yo = ^ ir + r")smtdt. 

cAo Jea 



'00 

Integrating by parts twice shows that 

xi9) -Xo = ^ {r{9fu'i9) - r(0o) V(^o)) 

CAo 

y{0) - Z/o = ^ {r{9fv'{9) - r{9orv'{9o)) , 

where 

u{9)=R{9) cos 9, v{9)=R{9)sm9 
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is the parametrization of the indicatrix used in Lemma 17^ Thus x{9) — Xq and y{9) — yo 
are simultaneously zero if and only if 



u'{9) 

v'{9) 



r{9y 



u'i9o) 
v'{9o) 



and because the indicatrix is strongly convex, this occurs precisely when 9 = 9o + 2mT for 
any integer n. Since 9 is unbounded as a function of arc length, it attains these values. 
Therefore, the projection of 7 onto the xy-plane is a simple closed curve. 

Finally, since dz = —^{xdy — ydx), Stokes' theorem implies that the difference z{9) — Zq 
along any horizontal curve in % is proportional to the signed area enclosed by the projection 
of the curve onto the xy-plane and the line segment connecting {x{9), y(^)) to (xq, yo). Thus 
z varies monotonically with increasing 6', and the projection of 7 onto the xy-plane is the 
curve of shortest Finsler arc length enclosing a given Euclidean area in the plane. D 

By the last observation in the proof above, the projections of sub-Finsler geodesies are 
precisely the solutions to the dual of the classical isoperimetric problem known as Dido's 
problem, named after Queen Dido in Virgil's Aeneid (see [MI). The classical solutions (using 
Riemannian arc length) are circles. In the Finsler case, the solution curves need not be 
circles, as one of the examples below illustrates. 

7.2.1. Randers metrics on "K and their geodesies. Consider the Randers-type, homogeneous 
sub-Finsler metric on D obtained by choosing the function r{9) to be 

r{9) = l + Bcos9, < B < 1. 

The indicatrix S^. for this metric in each plane D^ is the off-center ellipse with polar equation 

1 



R 



For this metric, 



l + 5cos^' 
3B sin 9 



2Vl + Bcos9' 
The geodesic equations ()7.4|) can be integrated explicitly in terms of 9. When Aq = 0, the 
geodesies are straight lines; when Aq 7^ 0, integrating yields: 



x{9) — Xq 



— (sin6' 
k 

1 



sin 6*0) 



(7.8) y{9) -yo = -{cos9~ cos 9q) 



z{9) - zo 



A(^ 



9o 



1 



sm 



9o)] + —[yo{sm9 



sin^n 



a;o(cos0 — cos 6*0)], 
These geodesies are liftings of circles of radius 1/k in the 



1 

2F 

where k = Xq a/(1 + Bcos9o) 
xy-plane. In the limiting case B = 0, we recover the geodesies of the fiat sub- Riemannian 
metric. When < i? < 1, the anisotropy of the indicatrix is manifested in the way that 
the area of the projected circle in the xy-plane (and, therefore, dz/d9) varies with the initial 
value 9o, unlike in the sub- Riemannian case. 
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Figure 1 shows some typical geodesies for 5 = |, starting from (xq, yo, Zq) = (0, 0, 0), with 
initial values Aq = 0.3 and ^o = 0, vr/2, and vr. For comparison, Figure 2 shows geodesies for 
the flat sub-Riemannian metric with these same initial values. 





Figure 1: Geodesies of the Randers 
metric 



Figure 2: Geodesies of the 
sub-Riemannian metric 



7.2.2. A "limagon metric" on 0-C and its geodesies. For an example in which the geodesies are 
not liftings of conic sections, consider the sub-Finsler metric whose indicatrix is the convex 
limagon with polar equation R = 3 + cos 6, so that 

r{9) ^ 



For this metric. 



3 + cos 9 



-3sin^(15cos^ + 13) 



2^/(9 cos^+ 11)3 
As always, the geodesies are straight lines when Aq = 0, but otherwise they are liftings of 
curves in the xy-plane defined by the equations 

1 /'sin6'(4cos6' + 6) sin6lo(4cos6'o + 6)' 
2L 



x{6) — xq 



(3 



cos( 



(3 



cos^n)^ 



y{e)-y, = -- 



9cos^+19 9cos^o 



19 



(3 + cos^)2 (3 + cos^o)^ 



where 



Ao 



V9cos^o + 11 



Ao^O. 



(3 + cos^o)" 

These curves in the xy-plane are not circles, nor are they ellipses (or limagons). Figure 3 
shows geodesies for this metric starting from (xq, |/o? -^o) = (0, 0, 0), with initial values Aq = 1 
and 9q = 0, 7r/2, and vr. Figure 4 shows the projections of these curves onto the xy-plane. 
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Figure 3: Geodesies of the limagon 
metric 



Figure 4: Projections of limagon 
metric geodesies onto the xy-plane 



8. Conclusion 

We have only begun to explore sub-Finsler geometry in this paper, and we have every rea- 
son to believe that it will become a useful extension of sub-Riemannian geometry, particularly 
in the context of control theory. In future papers, we plan to investigate higher-dimensional 
cases (including the important phenomenon of abnormal geodesies), singularities, and other 
topics likely to be of interest for control theory applications. 
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